Experimental targeting and control 
of spatiotemporal chaos in nonlinear optics 
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We demonstrate targeting and control over spatiotemporal chaos in an optical feedback loop 
experiment. Different stationary target patterns are stabilized in real-time by means of a two 
dimensional space extended perturbation field driven by an interfaced computer and applied in real- 
space to a liquid crystal display device inserted within a control optical loop. The flexibility of the 
system in switching between different target patterns is also demonstrated. 
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Control of complex dynamics refers to a process 
whereby the critical sensitivity of such dynamics to exter- 
nal disturbances is capitalized in order to select a proper 
tiny perturbation able to attain a desirable target behav- 
ior. In the last decade, a series of relevant issues such as 
stabilizing a given trajectory within an infinite set of un- 
stable periodic orbits embedded within a chaotic attrac- 
tor (control of chaos) or bringing a chaotic trajectory to 
a small neighborhood of some desired locations in phase 
space (targeting of chaos) have been addressed and solved 
in both low and high dimensional time-chaotic dynamics 



More recently, the interest switched to implement- 
ing control strategies for the stabilization of space-time 
chaotic dynamics occurring in extended systems. In this 
latter framework, some theoretical attempts and numer- 
ical demonstrations have been offered for achieving con- 
trol over one and two dimensional patterns Q, coupled 
map lattices and arrays of oscillators 0], or relevant 
model equations describing universal features of systems 
closeby to space time bifurcations, such as the Complex 
Ginzburg Landau Equation Q and the complex Swift- 
Hohenberg Equation yj. 

Robust and reliable experimental control over space 
time chaos (STC) remains however an open problem. In 
the field of nonlinear optics, a few experimental demon- 
strations of pattern control and targeting have been of- 
fered, based on Fourier filtering techniques 0,0- These 
methods provide a nice efficiency for the stabilization of 
stationary patterns with global symmetries, but their ap- 
plication to target patterns involving complex phase re- 
lationships among different Fourier components (such as 
localized patterns or selected snapshots of STC evolving 
patterns) is strongly limited by the practical difficulty of 
building suitable Fourier masks. 

In this Letter we show the first experimental evidence 
of control of STC based upon a real-space real-time feed- 
back technique, which is able to circumvent the above 
difficulties, thus allowing stabilization and targeting of 
two dimensional stationary patterns with arbitrary sym- 
metries and shapes. This is realized by means of per- 
turbing fields applied directly in the real-space and in 



times shorter than the characteristic time of the pattern 
dynamics, so that stationary non-homogeneous patterns 
of arbitrary complexity in space can be indifferently tar- 
geted and stabilized with good efficiency. We further- 
more show that our control strategy offers dynamical 
flexibility in switching from one to another target pat- 
tern, without the need of removing optical components 
(as e.g. filters) in the control loop. 

The experimental setup is sketched in Fig. la. It con- 
sists of a main optical feedback loop (MOFL) hosting a 
Liquid Crystal Light Valve (LCLV) [8|, and of an ad- 
ditional electrooptic control loop. The latter is essen- 
tially constituted by a videocamera, a personal computer 
driving a liquid crystal display (LCD), and a laser beam 
which traverses the LCD before being injected into the 
MOFL. The LCLV operates as a Kerr-like medium, i.e., 
it induces on the reading light a phase delay proportional 
to the writing intensity, over a wide range of input inten- 
sities. This proportionality relation holds for the experi- 
mental parameters used in the present investigation. The 
LCD display, operating in transmission, encodes linearly 
the gray level images output by the PC, onto the laser 
beam traversing it. 

When the control loop is open, a homogeneous wave 
is sent onto the front face of the LCLV, and is reflected 
acquiring a spatial phase modulation. The beam propa- 
gates within the MOFL from the front face of the LCLV 
to the input plane of a fiber bundle, experiencing diffrac- 
tion, and thus converting phase into amplitude modu- 
lations. The fiber bundle FB just relays the intensity 
distribution from the input to the output plane, which is 
in contact with the LCLV rear face. 

In these conditions, the dynamics within the MOFL is 
described by the equation |J E3 : 



dt 



4> ) + DV 2 (j) + aI fbl 



(1) 



where <fi(x, y) is the phase of the optical beam at the 
output of the valve, r is a relaxation time, D a diffusion 
coefficient, a the nonlinearity strength of the LCLV, <f>Q is 
the working reference phase, and Ifb(x, y) is the feedback 
intensity impinging at the rear side of the valve. 




(b) (c) 

FIG. 1: a) Experimental setup. Main loop: an extended laser 
beam is closed through a non-linear Kerr-like medium (liq- 
uid crystal optical valve). Instabilities develop in the trans- 
verse plane of the beam. M: mirror, Oi: microscope objec- 
tives; Pi: pinhole; A: aperture; BSi, BS2: beam splitters; 
LCLV: liquid crystal light valve; Li, L2: lenses of focal lens 
/; L5: additional lens; FB: fiber bundle. In our experiment, 
2/ — (h+fa) = +90 mm. Control arm: O2: microscope objec- 
tives; P2: pinhole; BS3: beam splitters; LCD: liquid crystal 
display; L3, L4: lenses. The arrows indicate the local direc- 
tion of light propagation, b) Snapshot of the uncontrolled 
STC state obtained for an input intensity I/I c ~ 3.2 (7 C be- 
ing the critical value for pattern formation from the uniform 
state to hexagons). The pattern intensity has been coded into 
a 256 levels gray scale, c) Space (vertical) -time (horizontal) 
dynamical evolution of the central vertical line of pixels of 
Fig. lb. 



The feedback intensity, which results from the propa- 
gation of the beam reflected by the front side of the valve 
trough the MOFL, is a nonlinear (and nonlocal) function 
of the phase <f> 0, El • On increase of the pump inten- 
sity J, the homogeneous solution destabilizes, resulting in 
hexagonal patterns close to threshold. If the pump is fur- 
ther increased, regular hexagons loose stability in favour 
of a space-time chaotic dynamics 0, |nj . Together with 
the pump value, another parameter of the utmost impor- 
tance is the spatial frequency bandwidth of the system 
Q, controlled by the aperture A in Fig. 1. Troughout 



(a) (b) (c) 

FIG. 2: Examples of target patterns (top row), controlled 
area in the system at 7 = 0.4 (center row) and corresponding 
far field images (bottom row) for the control trials of a perfect 
hexagonal pattern (a), a square pattern (b), and a snapshot 
of the uncontrolled dynamics (c) . 

the present paper, this bandwidth is kept fixed at ~ 1.3 
times the diffractive wavenumber of the system, which 
gives the scale of the unstable structures. 

In order to achieve control over the dynamics, a frac- 
tion of the beam traveling on the MOFL is extracted 
and detected by a videocamera, which is interfaced to 
the PC via a frame grabber. The computer processes the 
input image, and sends a suitable driving signal to the 
liquid crystal display device. The LCD transfer function 
T(x, y) can be written as the contribution of a constant 
mean transfer coefficient To, plus a modulation signal 
s(x,y). This modulation s is chosen to be proportional 
to the error signal between the current pattern intensity 
Ifb present in the system, and a target pattern Jt(x,j/): 

s(x, y, t) = -f3 (J/6 (a;, y, t) - I T (x, y)) . (2) 

The digital processing operations performed by the PC 
include the evaluation of the above error signal, and the 
calculation of the cross-correlation between pattern and 
target. The entire time series of pattern is also recorded 
on the hard disk for further off-line processing. 

The resulting refreshing time for the above procedure 
is at most 200 ms, to be compared with the characteristic 
time of the pattern dynamics (computed from the decay 
of the autocorrelation function) which is of the order of 
the second. 

The LCD is illuminated with a uniform intensity Jo, 
and the output beam is imaged onto the rear (writing) 
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FIG. 3: Correlation function C(t) (a,c,e) and amount of power 
p(t) injected within the control arm (b,d,f) vs. time (in sec). 
The two quantities are defined in the text. In all cases the 
continuous line refers to 7 = 0.4, the dotted line to 7 = 
0.2, and the dashed line to 7 = 0.1. (a,b) the target is a 
perfect hexagonal pattern; (c,d) the target is a snapshot of 
the uncontrolled dynamics; (e,f ) the target is a perfect square 
pattern. I/I c ~ 3.2. 



side of the valve. Following the above discussion, this 
beam will consist of a constant term T I that acts in 
renormalizing the valve working point (f>o to (f)Q , plus a 
modulated controlling beam sIq. 

The equation of motion when the control loop is closed 
is therefore: 



dt 



</)' ) + DA(l> + a(I fb -j(I fb -I T )). (3) 



where 7 = (3Iq. 

We initially set the light intensity at the input of the 
feedback loop to be I/I c ~ 3.2, where I c is the criti- 
cal value for pattern formation from the uniform state to 
hexagons. In these conditions, the uncontrolled evolution 
brings the system to display a time evolving, spatially 
disordered pattern, where many defects are continuously 
created and annihilated within an hexagonal-like pattern, 
thus generating STC . A typical snapshot of the un- 
controlled dynamics is reported in Fig. lb. Fig. lc 
reports the space-time dynamical evolution of the cen- 
tral vertical line of pixels in Fig. lb, showing how the 
uncontrolled dynamics evolves within STC, with a non 
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FIG. 4: Correlation function C(t) (see text for definition) 
vs. time during the sequential control trial at 7 = 0.4. The 
target patterns are a snapshot of the uncontrolled dynamics 
(10 < t < 30 sec), a square pattern (50 < t < 70 sec.) and a 
hexagonal pattern (90 < t < 110 sec). The dashed line indi- 
cates the switching on/off times. The vertical lines separates 
the three domains in time in which a different pattern is taken 
as target for the control. In each time domain, the correlation 
is calculated using the corresponding target pattern. 



stationary complex local dynamics and a decaying spa- 
tial correlation. 

Starting from these conditions, three different target 
patterns are selected, namely perfect hexagons, squares 
and a particular snapshot of the uncontrolled dynamical 
evolution (shown in the top row of Fig. 0). 

Perfect stationary hexagons are a stable solution close 
to the pattern formation threshold; they are destabi- 
lized when the pump is increased, resulting in the space- 
time chaotic dynamics here considered. Therefore, use of 
hexagons as a target tests the ability of the method to 
control an unstable solution. 

On the other hand, using as target a snapshot of 
the uncontrolled dynamics assesses the robustness of the 
method to freeze a given natural state of the uncontrolled 
dynamics, in the very same spirit as the so called target- 
ing of chaos . Finally, squares are never spontaneously 
selected by the system without control, and therefore 
they serve us to assess the ability of the control strat- 
egy to force the appearance of an arbitrary symmetry. 
The results obtained for 7 = 0.4 in the three cases are re- 
ported in the center row of Fig[21 indicating that the con- 
trol procedure is successful in all cases. The patterns are 
conveniently stabilized over a selected area of 128 x 128 
pixels (about 10 pattern wavelengths) in the system. The 
bottom row of Fig[2]shows the far field images of the con- 
trolled dynamics. While hexagonal and square patterns 
have a rather simple global symmetry (thus allowing for 
an easy implementation of Fourier filtering techniques, 
like the one of Ref. jfj), the target STC snapshot in- 
volves the presence of a complicate power spectrum. The 
fabrication of Fourier masks reproducing the amplitude 
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and phase of these patterns appears extremely difficult 
in experiments. 

Our control method circumvents such practical diffi- 
culties, and effectively stabilizes such complex Fourier 
patterns. Hence, it configures as the most reliable choice 
for the stabilization and targeting of two dimensional sta- 
tionary structures with arbitrary symmetries and shapes. 

To evaluate quantitatively the control ability of our 
method, we use the time dependent correlation function 
C(t) =< Ifb(r,t) ■ 2r(r) >r between the instantaneous 
pattern and the target one (< ... >r denotes a spatial 
average). We also measure the amount of power p(t) = 
j[< (Ifb(r,t) — Jr(r)) 2 >r] 1 / 2 injected within the control 
arm for controlling the target pattern. 

The correlations vs time are shown in Fig|3| for in- 
creasing values of 7 for the control task of a perfect 
hexagonal pattern (a), a snapshot of the uncontrolled 
dynamics (b), and a square pattern (c). At lower values 
of 7 (7 = 0.1,0.2), there is already a partial control of 
the dynamics, though several deviations of the patterns 
around the target ones still remain. This is reflected by 
the rather large fluctuations of C(t) visible in Fig. 3. The 
correlation value increases with 7 up to 7 = 0.4, when 
a high degree of control is achieved in all cases. As 7 
increases, the transient time before reaching control over 
the target pattern decreases. On the opposite, when con- 
trol is switched off, the relaxing time is determined by the 
characteristic time of the STC decorrelation, and is there- 
fore independent of 7. The value of the control power pit) 
needed to achieve control is in the range 0.1 — 0.2 for the 
data reported. 

Notice that the maximum correlation is attained for 
the case of the STC snapshot. This is due to the fact that 
in this case the target pattern is a specific configuration 
of the natural evolution of the dynamics. In the cases 
of squares and hexagons, instead, the target patterns are 
are digitally generated trying to reproduce at best the 
natural profile of the spots present in the system; but this 
procedure is intrinsically imperfect, because the exact 
profile of the spots is unknown. The mismatch between 
the profiles of the targets and the patterns lowers the 
correlation levels, with respect to the values that would 
be obtainable in the ideal case, in which the profile of the 
solution is perfectly reproduced by the target image. 

It is worth observing that the correlation when control 
is off in Fig. 3c does not decay to zero; this is due to 
the fact that our uncontrolled STC dynamics gives rise 
to a non zero mean field. This point can be qualitatively 
appreciated from inspection of Fig. lc: the pattern has 
a certain degree of "phase rigidity" , i.e., even if there are 
chaotic fluctuations, bright (dark) areas remain more or 
less bright (dark) for most of time. Similar properties 
have been observed experimentally and discussed in var- 
ious other cases of space extended systems giving rise to 
STC dynamics 

Finally, we demonstrate that our control technique is 



greatly flexible in dynamically switching between differ- 
ent target patterns. For this purpose, we prepare a time- 
dependent target pattern formed by the ordered sequence 
of a snapshot of the uncontrolled dynamics, a square pat- 
tern and a hexagonal pattern, each one presented to the 
system for a time T — 20 sec, and intermingled with 
periods of T = 20 sec. of time in which the system is 
left uncontrolled to reset the original state of STC. The 
results are shown in Fig. where one sees that the sys- 
tem is able to attain each one of the target pattern in 
the sequence for the same value of 7 = 0.4, as well as to 
switch between the different patterns. Notice that each 
target pattern in the sequence produces a different state 
in the correlation, depending upon its specific instabil- 
ity features within the uncontrolled STC regime. The 
maximal correlation is again obtained for the snapshot 
of the uncontrolled dynamics, since this pattern repre- 
sents a specific state compatible with the uncontrolled 
dynamics. 

In conclusion we have demonstrated targeting and con- 
trol over spatiotemporal chaos in an optical feedback loop 
experiment. Stationary target patterns with arbitrary 
symmetry and complexity are stabilized in real-time by 
means of a two dimensional space extended perturbation 
field at low cost in terms of the perturbation amplitude. 
The flexibility of the system in switching between differ- 
ent target patterns has also been demonstrated. 

Work partly supported by EU Contract HPRN- 
CT-2000-00158, and MIUR-FIRB project n. 
RBNE01CW3M-001. 
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